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Abstract 

The hadronization effects induced by different colour connections of ccqq system in e + e~ anni- 
hilation are investigated by a toy model where diquark fragmentation is employed based on Pythia. 
It is found that the correlations between the charm baryons and charm antibaryons produced via 
diquark pair fragmentation are much stronger, and their momentum spectra are harder than those 
from the standard colour connection in Pythia. 
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The hadronization for partons fragmenting into hadrons is very important, and can only be 
treated by hadronization models up to now. The colour connection plays a crucial role to set the 
"surface" between the Perturbative Quantum Chromodynamics (PQCD) phase and the hadroniza- 
tion one, while the topic related to the colour connection is beyond the approach of PQCD. It 
is argued that some clews can be drawn from analyzing the decomposition of the colour space 
of the final partons produced in PQCD phase [1-3]. For the colour space of the parton system, 
there are many decomposition ways corresponding to various colour connections, respectively 
(see, e.g., [4,5]). Hence, it is instructive to study the special properties of the final hadrons, which 
could give information for the evidence of specific colour connections. Such kinds of investigation 
have been made for various parton systems [6-8]. One interesting and important example is the 

<?2 <?3 <?4) system produced via hard collisions. This is the simplest system where many phe- 
nomena/mechanisms related to QCD properties could be studied, e.g., (re)combination of quarks 
in producing special hadrons, influence on reconstruction of intermediate particles (like V7 ± ) by 
soft interaction, etc., most of which more or less relate with the colour connections among the four 
(anti)quarks. In this letter, the possible colour connections of this system are briefly reviewed. It 
is pointed out that the hadronization effects of one kind of the colour connections have never been 
modeled yet. Its special case (diquark pair fragmentation) in e + e~ annihilation is investigated in 
detail. 

For qiqiqiq^ system, two of the decomposition ways of its colour space are given as follows: 

(3i®32)®(3 3 ®34) = (Ii2©8i2)®(l34©8 34 ) = (I12® I34) © (812® 834) e- ■ •, (1) 
(3i®3^)®(3 3 ®3 2 ) = (li4©8i 4 )®(l32© 8 3 2) = (li4<S)l32)© (814(8)832) ©■■•, (2) 
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Figure 1: (a) B — > n +D U and (b) B — > p ±E" processes. 

where 3, 3*, 1, 8 respectively denote the representations triplet, anti-triplet, singlet and octet of 
the Group SU C (3), and the subscripts correspond to the relevant (anti)quarks. In these two colour 
decompositions, quark and antiquark form a cluster/string which fragments into hadrons indepen- 
dently. This picture is adopted in the popular hadronization models (for detail, see [9-12]). The 
difference between these two decompositions leads to the colour reconnection phenomena as dis- 
cussed in the processes e + e~ — > W + W~ /zPzP — ^ qiqiq^qA - ► 4 jets at LEPII [2, 8, 13-16], etc. 
One might also notice, in the e + e~ — ► q\q\ + g* — > qxqxqiqi process, (q\q\) and (qiqi) cannot 
be in colour-singlets, while (#1^2) and (<?2<7i) are usually treated as colour- singlets and hadronize 
independently. This corresponds to the term 1 14 © I32 in eq. © [17]. 

The colour space of qiqiq^qA system can also be decomposed as follows: 

(3i © 3 3 ) ® (3* 2 © 3$) = (3j 3 © 613) © (3 24 © 6* 24 ) = (3f 3 © 3 24 ) © (613 © 6* 24 ) © ■ ■ • , (3) 

where 613 (624) denotes the representation sextet (anti-sextet) of the Group SU C (3). A special 
case is that two quarks in colour state 3* attract each other and form a "diquark" when their 
invariant mass is small enough, and similarly two antiquarks form an antidiquark. For this de- 
composition, only the whole four (anti)quarks system can form a colour singlet, hence they must 
hadronize altogether. This hadronization picture is quite different from that adopted in the pop- 
ular hadronization models. Unfortunately, this picture has not been taken into account though it 
is crucial and exists indeed in many processes. One example is the exclusive B decay as shown 
in fig. [T] where b — ► c + W~* — > c + d + u\. In this process, the quark-antiquark pairs du\ and cu 
are in colour-singlet, and naturally form two mesons and Do (figH^a)). Here the combination 
for du and cii\ is also possible. These two combinations correspond to the colour connections 
shown in eqs. and ©, respectively. The branching ratio of B — > n~D° has been measured 
to be (4.98 ±0.29) x 10~ 3 [18]. In the meantime, it is very interesting to notice that another 
kind of colour connection corresponding to eq. © exists, i.e., the cd and uu\ are in colour state 
3* and 3, and combine with a d or d from vacuum to form two baryons and p, respectively. 
This is a 'signal channel' for diquark- antidiquark pair fragmentation, and its branching ratio is 
(O.45+019 ±0.007 ±0.12) x 10~ 4 [19]. Additionally, the three body baryonic decay processes 
such as B -> A c pK- (Br = (2.1 ±0.7) x 10~ 4 [18]) can also be regarded to be induced by di- 
quark fragmentation. Another example is the doubly-heavy baryon production in high energy 
collisions [20,21]. Recently, SELEX Collaboration observed the hadronic production of S cc at 
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Figure 2: One of the four Feynman diagrams for e + e — > ccqq process. 

Fermilab [22]. This is dominated by the production of cc pair in PQCD process, and can also give 
an unambiguous confirmation of the existence of the diquark fragmentation. 

The aim of this letter is to investigate the diquark fragmentation in e + e~ — > qiqiqiqj — > h's 
process. The accumulated/accumulating data in various e + e~ colliders can provide sufficient 
statistics to study the hadronization of the four quark states. In order to specify different colour 
connections of the four quark system via measuring the corresponding hadron(s), we focus on the 
e + e~ — > ccqq — > h's process, where q represents u, d, or s quark. The ensemble properties related 
to a class of hadrons (e.g., charm baryons) are investigated, so that we can search for the quantities 
which could be measured in experiments. The statistics can also be enhanced by studying a group 
of hadrons. 

As mentioned above, the hadronization of the colour connection displayed in eq. © has not 
been considered in popular hadronization models. However, one can employ the event generator 
Pythia [10] based on Lund model [9], to describe the hadronization of the special case that both 
cq (in 3*) and cq (in 3) have small invariant masses so that they can be regarded as diquark and 
antidiquark, respectively. For this specific case, there are strong constrains on both colour state and 
phase space. Therefore it is important to investigate the fraction of this specific case, and relevant 
properties, at parton level. 

Parton level analysis At leading order QCD, the differential cross section for e+e —> ccqq 
process(fig. can be written as 

do = — d£ips4 \M\ 2 , (4) 
2s 

where s is the total energy square, L ips$ represents the 4-particle phase space and \M\ 2 is the spin 
averaged- summed invariant amplitude square. At the energy of B factories (eg., y/s = 10.52GeV), 
we calculate the joint distribution d 2 N/ [dm\dm'i) where m\ and ni2 respectively represent the in- 
variant mass of cq and cq. In our calculation, we take q = u as an example, and set the precise struc- 
ture constant to be 1/137, the strong coupling constant a s = 0.2, the quark masses m c = l.5GeV/c 
and m u c j = 0.33GeV /c 2 . The corresponding results are shown in fig. |3] It is important to notice 
that the production probability is not quite small when mi and m,2 are close to the diquark mass 5 . 
Especially, when the luminosity is large, the absolute event number with diquark pair produc- 
tion could be large enough to be measured. Employing the KEK B factory integrated luminosity 
(525.828/Z? 1 until Dec. 23, 2005) we obtain the number of the events with m\ . mi < 2.5 GeV /c 2 
to be around 3 x 10 4 . The number of events for q = d or q = s is at the same order. 

5 The diquark mass is around 2GeV /c 2 set by Pythia. 



Figure 3: The joint distribution d ^ ^ m for e + e — > ccww. 

It is well known that in the ^ y* — > // process, the polar angle distribution of the massless 
spin 1/2 fermion (J If) is predicted to be 1 + cos 2 B. For the 2-jet events initialized by e + e~ — > qq 
at energies much smaller than M z o, this will be approximately reproduced for the local parton 
hadron duality. For the two clusters cq and cq fragmentation, the event shape is very close to 2-jet 
event. If the polar angle distribution in this case is much different from the lowest order one, it can 
be expected as a possible signal to identify diquark pair fragmentation at B factory energy. The 
corresponding distributions with m\,mi< 2.5GeV are shown in fig. |4| Obviously the distribution 
is quite different from 1 + cos 2 Q, where refers to the angle between the cluster cq momentum 
and the e~ momentum. 

Before investigating the hadronization effects induced by different colour connections, we dis- 
cuss the theoretical uncertainties caused by quark mass, parton shower, etc. For charm quark mass, 
it is general to take the value between l.lGeV /c 2 to l.&GeV/c 2 , and for light quark mass, the 
same value as that in Pythia is adopted. For e + e~ — > ccqq (q = w, d) process, we calculate its cross 
section c(m q ) with respect to light quark mass m q . At the same time, we investigate the production 
cross section o(m c ) of e + e~ — > cccc. The ratios for a(m,)/a(m,o)(? = ^, c) withm^o = 0.33 GeV/c 2 
and m c Q = \ .5GeV /c 2 are shown in fig. |5](a). One can find that the cross section for e + e~ — > ccqq 
process is not quite sensitive to the light quark mass with charm mass fixed. This implies that 
the small light quark mass does not introduce any more dangerous singularities. The differential 
distributions l/ada/dM cq / cc are displayed in fig. |5](b). It is easy to find that the shape of the 
distributions for e + e~ — » ccqq and that for e + e~ — > cccc (fig. 13(b)) are quite similar, and a simi- 
lar study on doubly charm diquark fragmentation for four charm process might also be plausible. 
We also compare our results for the cross section of e + e~ — > cccc with those given in [23], and 
find completely agreement. The cross section for four quark production in e + e annihilation at 
high energies (eg., LEP1 and LEP2) may be affected by the gluon radiation, which is simulated 
by parton shower process in Pythia. Fortunately, at B factory energies, y/s ~ \0GeV, the largest 
energy of each quark is about 5GeV. If it could radiate gluons further, its virtuality will quickly 




Figure 4: Normalized polar angle distribution of diquark cq (solid histogram). The dashed curve 
is for the normalized 1 + cos 2 Q distribution. 



reach the lower limit that perturbative QCD can be applied safely. As a result, the parton shower 
process at B factory energies would not play an important role as that at LEP energies. When one 
investigates physics for the e + e~ — > h's process at LEP or even higher energies, the parton shower 
process must be taken into account. Unfortunately, till now, it is still unclear how to describe the 
gluon radiation from diquark pair in annihilation. As the first step of the investigation, we 
concentrate on the physics at the B factory energy, where the influence of gluon radiation seems to 
be small and could be neglected. 

Hadronization For the — > ccqq process, when m\ and mi near the mass threshold 
of diquark, we employ Pythia to investigate the differences in final hadron states corresponding 
to two kinds of colour connections. One is that cq and cq form diquark-antidiquark pair, and 
hadronize altogether. In the picture of Lund model, diquark and antidiquark are in 3* and 3 colour 
state, respectively. As a result, one colour string will be stretched between them, and the string 
will fragment into hadrons at the end. The other case is that cq and cq form two colour singlet 
clusters (strings). These two strings will fragment into hadrons independently. This hadronization 
picture is widely adopted in the popular event generators, and hereafter we denote it as '2- string 
fragmentation'. The hadronization of these two cases can be treated by the corresponding subrou- 
tines in Pythia. For simplicity, we fix the invariant masses of the (anti-)quark pair cq (cq) to be the 
(anti)diquark mass used in Pythia. The results obtained under this simplification can be considered 
as a good approximation of the average over a reasonable (hence small) range of cqlcq cluster 
mass. It is straightforward to adopt that the quantities we calculate in this section is not sensitive 
to the cluster mass in a small range around the diquark mass set by Pythia. By a similar argument, 
we also assume the inner orbital angular momentum of the cq (cq) cluster to be 0. However, for 
treating the second kind of colour connection, the inner relative movements between c (c) and q 
(q) are still not fixed and generated with the relative weight given by the differential cross section 
at parton level discussed above. 





Mcq(c) (GeV/O 



Figure 5: (a) The dependence of cross section on quark mass for e + e~ — > ccqq process (solid 
line) and e+e~ — > cccc process (dashed line), (b) The normalized distribution 1/a dc/dM for 
e + e~ — > ccqq process (solid line) and e + e~ — > cccc process (dashed line). 





Figure 6: (a) X p (= with p max = J f—m 2 ) spectra of A c at ^fs = \0.52GeV, from diquark 

fragmentation (solid line), 2-string fragmentation (dashed line) and e + e~ —> cc — > A c +X (dotted 
line), (b) Comparing with data [24]. The histograms correspond to three different baryon produc- 
tion scenarios in Pythia: solid line for 'diquark', dotted line for 'simple popcorn' , dashed line for 
'advanced popcorn'. 



In the diquark pair fragmentation process, it is straightforward that diquark tends to fragment 
into baryon, which is in general the leading particle. If the charm baryons from the diquark frag- 
mentation are detected, their momentum spectra should be harder than those from the 2-string 
fragmentation, or those from the general — > cc — > h's process. This feature can be displayed 
by the momentum (fraction) spectra of A f at y/s = 10.52GeV (fig. El(a)). Recently, the momentum 
spectra of various charm hadrons have been measured by BELLE Collaboration [24]. To make full 
use of the available data and to identify whether the diquark fragmentation events exist or not, we 
must make clear to what extent Pythia can reproduce the data. For charm meson production, we 
find that the predictions of Pythia agree to the available data quite well (also see [25]). For our 
aim, we concentrate on the investigation of charm baryon (eg., A c ) production. In Pythia, there 
are three scenarios to describe baryon production, i.e., diquark, simple popcorn and advanced pop- 
corn. Their predictions can be consistent with each other provided that the relevant parameters are 
tuned (fig. E](b), the parameters all take default values except PARJ(\) = 0.2, PARJ(\S) = 0.19 
and PARF{\92) = for advanced popcorn). From fig. |6](b), one can find that the distributions of 
A c predicted by standard Pythia are slightly softer than the data, while the momentum spectrum 
obtained by diquark fragmentation is slightly harder. It is also interesting to notice that for the di- 
quark fragmentation process, the small momentum region is slightly enhanced for the mass effect 
(fig- ED(a)). This phenomenon is also indicated by the data shown in fig. 0b), though the statistical 
error is too large to give a definite conclusion. Because of the phase space limit, the probability 
of diquark pair production may not be large. As a result, the predictions for A c production includ- 
ing the diquark pair production will not conflict with the available data, and a portion of diquark 
fragmentation events is favoured. 

For further investigation, we calculate thrust distribution 1 /N dN/dT for the diquark fragmen- 
tation, 2-string fragmentation and e + e~ — > cc — > h's process (fig. Eta)). One can find that the thrust 
distribution of the diquark fragmentation is much 'thinner' than those of the other cases (fig. Wia)). 
Aware that the diquark fragmentation events are mostly 2-jet like, we compare the thrust distribu- 
tions of 2-jet events for the three cases. We use the JADE algorithm with y cut = 0.12. In this case, 
around 90% diquark fragmentation events are 2-jet. The properties are similar to those without 
jet constrain (fig. 13b)). According to these discussions, further experimental investigation of the 
colour connections shown in eq. © can be done by selecting 2-jet charm events at the B factory, 
and choosing those with large thrust (say, larger than 0.85). 

It is natural that in the fragmentation of diquark-antidiquark pair, a baryon is produced together 
with an antibaryon. That means, the correlation between charm and anti-charm baryons in the 
diquark pair fragmentation is much stronger than that in the single charm quark fragmentation 
(either the 2-string fragmentation or the general e + e~ — > cc — > h's). We know that, in the latter two 
cases, the production probability for charm meson is much larger than that for charm baryon. This 
phenomenon can be illustrated by the baryon (B) antibaryon (B) correlation, defined as, 

, _. 2N(BB) 
A(BB) = — — y - — (5) 
V ; N(B)+N(B)' 

where N(BB), Nb and N§ respectively denote the probability for finding BB pair, B and B. We 
take the correlation of A c and A c as an example in Table [l] One can notice that the correlation 
in the diquark pair fragmentation process is stronger than the other cases by about 10 times. This 
correlation is not sensitive to the gluon radiation of the diquark pair so we also give the results at 
Z° pole. This quantity could be a good probe for the diquark pair fragmentation. The results for 
other charm baryons, such as E c , S c and the corresponding spin-3/2 particles are similar. 

To summarize, we investigate the hadronization effects induced by different colour connections 
of four quark system (eg., ccqq) in e + e~ annihilation. It is interesting to point out that besides the 
normal colour structure with two colour singlets cq and qc which fragment into hadrons indepen- 
dently in the DODular models, there exists another kind of colour structure, as shown in ea. (l3l). 




Figure 7: Thrust distributions of three cases at y/s = \0.52GeV . (a) Diquark fragmentation (solid 
line), 2-string fragmentation (dashed line) and — > cc — > /j's (dotted line), (b) Same as (a) 
except that we select 2-jet events, with JADE Algorithm and y cut = 0. 12. 





e + e — > cc — > /I'i 


Diquark Fragmentation 


2-string 


10.52GeV 


0.057 


0.76 


0.044 


91.19GeV 


0.081 


0.72 


0.081 



Table 1: A C A C flavor correlations. 



which has to be taken into account in some cases. For that colour connection, needless to say, the 
hadronization of these four quarks should be treated as a whole system, i.e., they interact among 
each other during the hadronization process. In this letter, considering an extreme limit case that 
cq and cq form diquark-antidiquark pair, we propose a toy model based on Pythia to describe the 
hadronization of them. It is found that the results originated from diquark fragmentation are sig- 
nificantly different from those of the ordinary colour connections. However because of the phase 
space limit, the hadronization effect from diquark pair fragmentation is not very large. To un- 
derstand the hadronization mechanism induced by different colour connections of final partons 
produced perturbatively, further theoretical and experimental investigations are still necessary. 
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